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ABSTRACT 

I present the results of a programme to monitor the strengths of the Li I 6708A, 
Kl 7699A and chromospheric Ha lines in a group of cool Pleiades stars. Consistent 
instrumentation and analysis techniques are used to show that there is no Li I vari- 
ability on timescales of 1 year that could possibly account for the apparent spread in 
Li abundances seen in Pleiades stars between effective temperatures of 4800-5200 K. 
Comparison with published data reveals tentative evidence for variability on 10 year 
timescales, but at a very low level. The lack of chromospheric activity variability above 
levels of 20 to 30 percent makes it difficult however, to rule out evenly distributed mag- 
netic activity regions causing a scatter in the Li I line strengths at a given abundance. 
The similar star to star scatter of K I line strengths in these and published data rein- 
forces the conclusion that it is still unsafe to attribute the Li I line strength dispersion 
to a large variation in Li depletion at a given mass. 
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1 INTRODUCTION 

Lithium is unique among metals - it is produced in the big 
bang, but is destroyed by p,a reactions at low temperatures 
{^ 2.4 X 10^ K) in stellar interiors. Observing Li in the cool 
stars of young open clusters is an excellent way of explor- 
ing and calibrating the physics of internal stellar mixing due 
to convection or additional "non-standard" processes. Rapid 
angular momentum loss as stars reach the ZAMS may make 
Li depletion due to poorly understood mechanisms such as 
rotational mixing, much more dramatic. Such investigations 
are important in understanding stellar physics, but also in 
deciding whether present day Population II Li abundances 
have hardly changed from the primordial big bang nucle- 
osynthesis value. This is supported by "standard" models 
that feature only convective mixing (Deliyannis, Demarque 
& Kawaler 1990, Bonifacio & Molaro 1997), but other mod- 
els including "non-standard mixing" predict that the Popu- 
lation II Li has been depleted from primordial values a factor 
of 2-3 higher (Deliyannis & Ryan 1997, Pinsonneault et al. 
1998). 

A particular problem for the standard models that has 
emerged in the last 15 years is the observation that cool 
stars (4500 K< Te// <5400 K) in the Pleiades (age 100 Myr) 
and Alpha Per (age 50 Myr) clusters exhibit an apparent Li 
abundance spread of more than an order of magnitude at 
a given mass (Duncan & Jones 1983, Butler et al. 1987, 
Soderblom et al. 1993a [hereafter - S93a], Balachandran, 



Lambert & Stauffer 1996, Randich et al. 1998). This spread 
cannot be reconciled with standard models, which predict a 
unique amount of pre-main sequence (PMS) Li depletion at 
a given mass and age. 

Several explanations have been put forward which as- 
sume that the observed spread corresponds to a real disper- 
sion in the Li abundances. These have been guided by the 
correlation (albeit an imperfect one) of reduced Li depletion 
and rapid rotation (S93a, Martin 1997, Randich et al. 1998). 
Extra mixing as stars spin down towards the ZAMS is one 
possibility (Chaboyer et al. 1995). The observed spread in 
rotation rates at the ZAMS could lead to different efficien- 
cies for rotationally induced mixing. Unfortunately, even the 
most recent models, which incorporate early angular mo- 
mentum loss via a circumstellar disk, seem unable to ex- 
plain the size of the abundance dispersion (Pinsonneault et 
al. 1998). Others have postulated that structural changes as- 
sociated with rapid rotation (Martin & Claret 1996) or dy- 
namo induced magnetic fields in the convection zone (Ven- 
tura et al. 1998) might reduce the efficiency of standard 
PMS Li depletion in rapid rotators. Other combinations of 
mechanisms have been plausibly suggested, such as mixing 
by internal gravity waves plus the magnetic blocking of me- 
chanical flux through the bottom of the convective zone in 
rapid rotators (Schatzman 1993, Montalban & Schatzman 
1996). These ideas seem promising because the observed Li 
depletion is significantly less than the conservative predic- 
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tions due to standard PMS convective mixing alone (Jeffries 
& James 1999). 

Another school ol thought is that the apparent Li abun- 
dance dispersion might largely be explained in terms ol a 
scatter in the equivalent width (EW) ol the Li I 6708A line 
(almost exclusively used lor this work) at a given Li abun- 
dance - caused by NLTE, chromospheres, or large-scale pho- 
tospheric inhomogeneities such as plages or starspots. It has 
long been known that the reduced temperature in sunspots 
increases the observed EW ol the Li I 6708A hne (Giampapa 
1984). Young, cool stars are known to have much larger 
photospheric inhomogeneities than the Sun, which are as- 
sociated with their enormous magnetic activity. The most 
recent modelling work (Carlsson et al. 1994, Houdebine & 
Doyle 1995, Stuik, Bruls & Rutten 1997) indicates that these 
effects may indeed induce a scatter in the apparent Li abun- 
dances il, as is normally the case, homogeneous, one dimen- 
sional atmospheres and Tefi-colour relationships are used to 
convert observed EWs into abundances. Impetus has been 
given to these explanations by two pieces of observational 
work. S93a (see also Stuik et al. 1997) find that the Kl 
7699A line, formed in similar conditions to the Li I line, also 
shows a rotation-dependent EW scatter (albeit smaller than 
that exhibited by Li I) at a given colour. As no potassium 
abundance variations are expected in the Pleiades, this re- 
sult challenges our understanding of how neutral alkali lines 
are formed in cool stellar atmospheres. Russell (1996) ex- 
amined the subordinate line of Li I at 6104A in several late- 
type Pleiads. He argued that Li abundances from this line 
are lower and show considerably less scatter than from the 
Li I 6708A line, possibly because of different sensitivities to 
NLTE conditions and chromospheric activity. Russell's con- 
clusions have since been challenged because of the difficulty 
in deblending the weak 6104A line from stronger nearby lines 
(Martin 1997). 

A useful test of whether Li abundance variations are 
real is to search for variability in the Li I 6708A EW. If the 
scatter in Li abundances is caused by large-scale atmospheric 
inhomogeneities one might expect to see EW variations as- 
sociated either with magnetic activity cycles or with the 
rotational modulation of spots and plages. Spots are often 
demonstrably asymmetric in their surface distribution on 
young, active stars, which for instance, allows rotation pe- 
riods to be measured from light curve modulations (O'Dell 
& Collier-Cameron 1993). The search for Li EW variabil- 
ity has a chequered history and is reviewed by Fekel (1996). 
There are individual examples of young active field stars 
where EW variations on timescales of the rotation period 
have been reported {e.g. Robinson, Thompson & Innis 1986, 
Basri, Martin & Bertout 1991, Batterer et al. 1993, Jeffries 
et al. 1994, Elbe et al. 1999) and other studies which claim 
to detect no variations at the 5% level or less {e.g. Boesgaard 
1991, Pallavicini et al. 1993). The only study carried out so 
far in young clusters is that of S93a. They say that obser- 
vations of some cool Pleiades stars taken on several occa- 
sions show no EW variations beyond their estimated errors 
of ~20mA, but no details are given. As the Li I 6708A EWs 
in late-G and early-K type Pleiads range from ~ 100 mA to 
~300 mA this amount of variability would seem unable to 
account for any significant fraction of the apparent Li abun- 
dance dispersion. 

To put this important result on a firmer footing, I have 



begun a programme to monitor the Li I 6708A line in a selec- 
tion of late-G and early-K type Pleiads over long timescales; 
to test for variability, to see whether any variability is also 
present in the K I 7699 A line and to see whether any variabil- 
ity is correlated with changes in chromospheric activity. In 
this short paper I present the first results of this programme 
which allow a precise investigation of variability on 1 year 
timescales using data taken with identical instrumentation 
and on ~10 year timescales for the first time, by comparison 
with the results of S93a. 



2 OBSERVATIONS AND ANALYSIS 



2.1 Target selection 

A dozen Pleiades stars were selected for re-observation from 
Table 1 in S93a. They were chosen to have 4800 <Tcff < 
5200 K, based upon the colour- Tcff relation used by S93a 
in their analysis, were not known spectroscopic binaries and 
had a large range of projected equatorial velocities (we sin i) 
and Li I 6708A EWs. 



2.2 Spectroscopy 

The spectroscopy was performed at the 2.5-m Isaac New- 
ton Telescope (INT) on the nights of 1997 November 10 and 
1998 November 26. In both cases precisely the same instru- 
mental set up was used. The Intermediate Dispersion Spec- 
trograph, 500 mm camera, H1800V grating and a TEK 1024 
pixel square detector were used to cover the spectral region 
from 6510A to 6755A at a dispersion of 0.23A per 24 pm 
pixel. A 1 arcsec slit projected to just over 2 pixels on the 
CCD, giving an instrumental spectral resolution (confirmed 
by arc lines) of 0.48A. A second wavelength setting was also 
used that covered the range 7570A to 7790A at 0.22A per 
pixel and an instrumental resolution of 0.46A. 

The usual calibration bias frames, tungsten fiat fields 
and B star spectra were taken, as well as copper/neon 
lamp exposures at every target position. A set of "mini- 
mum chromospheric activity" standard stars, from the list 
in Soderblom et al. (1993b), were also observed at high signal 
to noise (S/N). The data were reduced and spectra extracted 
and wavelength calibrated using the Starlink FIGARO soft- 
ware package. Two exposures of each Pleiades target were 
taken at each wavelength setting. The exposure times were 
judged to yield S/N levels of 80-110 per pixel in the contin- 
uum of the co-added Li I 6708A spectra and 55-70 per pixel 
in the continuum of the co-added Kl 7699A spectra. 

The seeing was not particularly good on either night 
(1.5-2.5 arcsec) and the final 4 hours of the 1997 observa- 
tions were badly cloud affected. Nevertheless I was able to 
complete observations of six of the programme targets in 
1997 and nine in 1998. A log of the observations and the 
S/N levels achieved in each co-added spectrum is given in 
Table 1. 

EWs of the Li I 6708A, Cal 6718A and Kl 7699A lines 
were measured by direct integration below continuum lev- 
els which were determined from polynomial fits to approx- 
imately line free regions of the spectra. For targets where 
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Table 1. A log of the INT spectroscopic observations and EW measurements. 



Name ( 


B-V)o 


Ve sin i 


Date 


UT 




AA6510 - 6755 




AA7570 - 7790 












S/N 


LiI6708A 


CaI6718A 


Excess Ha 


S/N 


K I7699A 






kms~^ 








EW (mA) 


EW (mA) 


EW (mA) 




EW (mA) 


Hz 34 


0.89 


< 7 


10 Nov 1997 


20:30 


90 


154 ±7 


186 ±6 


158 ± 12 


60 


261 ± 13 


Hz 345 


0.81 


18 


10 Nov 1997 


21:00 


80 


282 ± 10 


176 ±7 


1508 ± 15 


60 


336 ± 12 


Hz 174 


0.81 


28 


10 Nov 1997 


21:30 


90 


299 ± 12 


220 ± 14 


1770 ± 30 


60 


391 ± 24 


Hz 916 


0.83 


< 7 


10 Nov 1997 


22:15 


95 


187 ±9 


180 ±7 


217 ±15 


60 


282 ± 11 


Hz 2407 


0.91 


< 7 


10 Nov 1997 


23:00 


100 


116 ±7 


218 ±6 


340 ±9 


65 


317 ±12 


Hz 2034 


0.93 


75 


11 Nov 1997 


00:00 


85 


220 ± 13 


238 ± 14 


1060 ± 50 


70 


469 ± 16 


Hz 345 


0.81 


18 


26 Nov 1998 


20:30 


100 


270 ±5 


172 ±6 


1250 ± 15 


60 


337 ± 15 


Hz 174 


0.81 


28 


26 Nov 1998 


21:00 


90 


290 ± 11 


203 ± 14 


1750 ± 50 


55 


433 ± 32 


Hz 916 


0.83 


<7 


26 Nov 1998 


22:00 


110 


185 ±4 


169 ±6 


207 ± 12 


60 


261 ± 15 


Hz 2407 


0.91 


<7 


26 Nov 1998 


23:00 


95 


125 ±5 


210 ±8 


351 ±9 


55 


324 ± 16 


Hz 2034 


0.93 


75 


27 Nov 1998 


00:00 


80 


230 ± 19 


253 ± 13 


1490 ± 80 


60 


452 ± 20 


Hz 34 


0.89 


< 7 


27 Nov 1998 


01:30 


110 


146 ±5 


191 ±5 


169 ± 10 


60 


277 ± 10 


Hz 263 


0.84 


10 


27 Nov 1998 


02:40 


100 


240 ±5 


166 ±6 


359 ± 11 


65 


300 ± 10 


Hz 1095 


0.86 


<7 


27 Nov 1998 


03:30 


90 


130 ±6 


187 ±9 


180 ± 13 


55 


285 ± 15 


Hz 1124 


0.91 


7.5 


27 Nov 1998 


04:30 


90 


189 ±7 


253 ±9 


527 ± 15 


60 


335 ± 12 



I have a pair of observations from 1997 and 1998, I was 
careful to use identical continuum regions, after correcting 
the spectra to a similar rest velocity, so that EW measure- 
ments are directly comparable for these stars. Note that it 
was not possible to select identical continuum regions for all 
the targets because the rapid rotators suffer from extensive 
blending and consequent blanketing of the continuum. The 
K I line is affected by a small amount of telluric absorption. 
This was corrected for by reference to a B-star spectrum. 

I also measured the Hq feature as an estimator of chro- 
mospheric activity. This is usually in absorption in inac- 
tive stars, but can be filled in or even show emission in the 
most active G/K stars. Again, consistent continuum mea- 
surements were used in pairs of observations to aid compar- 
ison. Following Soderblom et al. (1993b), I subtracted the 
spectrum of each Pleiades star from the spectrum of a nor- 
malised minimum activity star of similar intrinsic colour. 
The standards used are the same as those in Soderblom 
et al. (1993b), namely, HD 3651, HD 4256, HD 4628, HD 
10476, HD 166620 and HD 182488. For the two stars with 
Ve sin i> 20kms^^, the standard star was broadened to the 
appropriate Ve sin i, before subtraction. The excess EW of 
the resultant "emission line" can be taken as a measure of 
the chromospheric activity and is comparable with the mea- 
surements in Soderblom et al. (1993b). Again, telluric fea- 
tures were visible in some of the spectra and were removed 
by reference to a scaled high S/N spectrum of a rapidly ro- 
tating B-star. 



3 RESULTS 

The EW measurements for all the relevant lines are given 
in Table 1, along with the S/N in the co-added spectra. 
Tables 2-5 list these results again and compare them with 
the results of S93a and Soderblom et al. (1993b). In Ta- 



Table 2. A comparison of deblended Li I 6708A EWs from this 
paper and Soderblom et al. (1993a). 



Name Deblended Li I EW (mA) 

S93a Nov 1997 Nov 1998 



Hz 34 


134 ± 18 


139 ±7 


131 ±5 


Hz 174 


241 ± 25 


286 ± 12 


277 ±11 


Hz 263 


290 ± 12 


- 


226 ±5 


Hz 345 


245 ± 12 


269 ± 10 


257 ±5 


Hz 916 


199 ± 12 


173 ±9 


171 ±4 


Hz 1095 


138 ± 12 


- 


116 ±6 


Hz 1124 


217 ±18 


- 


174 ±7 


Hz 2034 


222 ± 36 


204 ± 13 


214 ± 19 


Hz 2407 


125 ± 12 


101 ±7 


110 ±5 



Table 3. A comparison of Cal 6718A EWs from this paper and 
Soderblom et al. (1993a). 



Name 




Cal EW (mA) 






S93a 


Nov 1997 


Nov 1998 


Hz 34 


186 ± 18 


186 ±8 


191 ±5 


Hz 174 


269 ± 25 


220 ± 14 


203 ± 14 


Hz 263 


218 ±12 


- 


166 ±6 


Hz 345 


182 ± 12 


176 ±7 


172 ±6 


Hz 916 


203 ± 12 


180 ±7 


169 ±6 


Hz 1095 


216 ±12 


- 


187 ±9 


Hz 1124 


237 ± 18 


- 


253 ±9 


Hz 2034 


256 ± 38 


238 ± 14 


253 ± 13 


Hz 2407 


229 ± 12 


218 ±6 


210 ±8 
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Table 4. A comparison of K I 7699A EWs from this paper and 
Sodcrblom ct al. (1993a). 



Name 




K I EW (mA) 






S93a 


Nov 1997 


Nov 1998 


Hz 34 


233 ± 18 


261 ± 13 


277 ± 10 


Hz 174 


- 


391 ± 24 


433 ± 32 


Hz 263 


250 ±12 


- 


300 ± 10 


Hz 345 


321 ±12 


336 ± 12 


337 ± 15 


Hz 916 


222 ± 12 


282 ± 11 


261 ± 15 


Hz 1095 


235 ± 12 


- 


285 ± 15 


Hz 1124 


369 ± 18 


- 


335 ± 12 


Hz 2034 


- 


469 ± 16 


452 ± 20 


Hz 2407 


294 ± 12 


317 ±12 


324 ± 16 



Table 5. A comparison of excess Ho EWs from this paper and 
Soderblom et al. (1993b). Where several measurements were given 
by Sodcrblom ct al. (1993b), all are listed. 



Name 


Excess Ha EW (mA) 




S93b 


Nov 1997 


Nov 1998 


Hz 34 


230 ± 15 


158 ± 12 


169 ± 10 


Hz 174 


- 


1770 ± 30 


1750 ± 50 


Hz 263 


480 ± 15 
340 ± 15 
390 ± 15 




359 ± 11 


Hz 345 


1570 ± 15 
1050 ± 15 
1240 ± 15 


1508 ± 15 


1250 ± 15 


Hz 916 


280 ± 15 
290 ± 15 


217 ±15 


207 ±12 


Hz 1095 


140 ±15 
340 ± 15 


— 


180 ± 13 


Hz 1124 


580 ± 15 


- 


527 ±15 


Hz 2034 


- 


1060 ± 50 


1490 ± 80 


Hz 2407 


230 ± 15 


340 ±9 


351 ±9 



ble 2 I compare the "deblended" Li I 6708A EWs quoted by 
S93a. The blend in question is with a weak Fel feature at 
6707.44A, which has a strength of 20(B - V)o - 3 mA ac- 
cording to S93a. The EWs listed in Table 2 have all been 
corrected according to this formula. Figures 1 and 2 show 
examples of the spectra from 1997 and 1998 for compari- 
son of the Li I, Ca I and K I lines. Figure 3 illustrates the 
spectral subtraction technique used for Ha and shows an 
example of the comparison between excess Ha emission in 
1997 and 1998. 



3.1 Variability on one year timescales 

The strength of the small dataset presented here is that 
the observations were taken with identical instrumentation 
and reduced in a very consistent manner. Comparison of 
the six Pleiades stars in common to the November 1997 and 
November 1998 datasets allows a strong statement on the 
amount of variability on yearly timescales, independent of 
any possible systematic errors in the EW measurements. It 
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Figure 1. Spectra for Hz 174 and Hz 916 around the Li I 6708A 
line from 1997 (dots) and 1998 (solid lines). 
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Figure 2. Spectra for Hz 916 and Hz 2034 around the KI 7699A 
hnc from 1997 (dots) and 1998 (solid lines). 



is clear from Tables 2-4 that the EWs of the Li I Ca I and K I 
lines are completely in agreement. For Li I the biggest (and 
most significant) discrepancy is (12±11) mA for Hz 345. The 
weighted mean discrepancy for all six stars is (2 ± 4) mA, 
with an rms discrepancy of less than 14 mA at a 95 percent 
confidence level. 

The EW discrepancies can be somewhat larger for the 
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6550 6560 6570 

Wavelength (A) 



6580 



Figure 3. The bottom panel shows the spectral subtraction tech- 
nique applied to Hz 345 (dots) around the Hq line in 1998. The 
template star (solid line) is HD 182488. The top panel shows the 
subtracted "excess" Ha emission lines for Hz 345 in 1997 (dots) 
and 1998 (solid line). 



continuum placement may be to blame for most of the ap- 
parent variation. From Tables 2 and 3, we can see that this 
may be the case for Hz 263, Hz 916, Hz 1095 and Hz 2407, 
but Hz 174 and Hz 1124 show Cal EW discrepancies in the 
opposite sense to those in Li I. If all the statistical error es- 
timates are correct, then we have to conclude that there is 
some genuine Li I EW variation at the level of ~ 30 mA in 
at least some Pleiades stars on ~ 10 year timescales. 

The evidence for variation in the Kl line is also diffi- 
cult to find. There is some indication that the S93a EWs 
are systematically lower than those presented here, which 
make the results for Hz 1124 intriguing. It seems that the 
EW was larger in the S93a observation, perhaps behaving 
in a similar fashion to the Li I EW. The other six stars in 
common, would be quite comparable if I have generally over- 
estimated the continuum level around the Kl line by ~ 3% 
with respect to S93a. Differences of this order are just plau- 
sible because of the many telluric features close to the line, 
the different resolutions of the observations and perhaps an 
under-estimation of scattered light in S93a's echelle obser- 
vations. 

Lastly, we can look at chromospheric activity. 
Soderblom et al. (1993b) list a number of measurements for 
some stars taken in 1988 and 1990, and our additional mea- 
surements paint a consistent picture of low-level variability 
of less than a factor of two on timescales of 1-10 years. In- 
terestingly, Hz 1124 is the most constant star in our sample 
(in terms of percentage variation) in contrast to the results 
from the other lines. 



other two lines, although their errors are larger and the sig- 
nificance of the discrepancy is still less than 1.5a in all cases. 
From the excess Ha EWs, there is clear evidence for vari- 
ability in the chromospheric emission of Hz 345 and Hz 2034 
at the 20-30 percent level, whereas the other four stars are 
constant to better than 10 percent. 



3.2 Variability on ten year timescales 

The dataset provided by S93a allows variability on longer 
timescales to be probed. Here we have a slightly larger 
dataset (nine stars), but there could be systematic differ- 
ences in the EW measurements due to choice of continuum 
placement and slightly different spectral resolutions. Ignor- 
ing these possible systematic errors for the present, we can 
see that the 1997/98 Li I line measurements are remarkably 
consistent with those of S93a taken in 1988 and 1990. The 
maximum and most significant discrepancy is (64 ± 13) mA 
for Hz 263. The weighted mean discrepancy between S93a 
and our 1998 observations is (21 ± 5) mA, in the sense that 
the S93a EWs are larger. However, from the error estimates, 
this cannot be explained (is statistically unacceptable) as a 
systematic continuum placement problem and there may be 
genuine variation in the Li I EWs at an rms level of about 
30 mA. 

A useful check on any systematic error is provided by 
the nearby Ca I line, which will have been measured with re- 
spect to a similar adjacent continuum region. The Cal line 
is formed from an excited level and has been shown to be 
considerably less sensitive to photospheric inhomogeneities 
than the Li I line. If both the Li I and Cal line show EW 
discrepancies in the same sense and of similar size then the 



4 DISCUSSION 

The new measurements of Li I EW presented here were ob- 
tained with identical instrumentation and taking care to 
avoid any systematic errors. The conclusion drawn from a 
small sample is that the Li I EW is not variable on one year 
timescales at more than a level of about 10-15 mA. For a 
late G/early K star with an effective temperature of 5000 K 
the curves of growth presented in S93a can be used to esti- 
mate that this would produce a spread in Li abundance of 
< 0.1 dex for the most Li abundant stars in the Pleiades and 
perhaps < 0.15 dex for the least Li abundant Pleiades stars 
at this temperature. Variations at this level cannot there- 
fore explain the vast majority of the more than one order 
of magnitude apparent Li abundance scatter in the Pleiades 
late G/early K stars. 

On longer timescales of ~ 10 years, the data do sup- 
port a limited amount of Li I variability, perhaps as much 
as an rms of 30 mA. While this is still far too small (0.2- 
0.3 dex) to explain the Li abundance scatter, it does agree 
with the kind of Li I variations which have been seen on 
some, but not all, young, nearby, rapidly rotating K-type 
Pleiades analogues. For example Hussain et al. (1997) mea- 
sure a 45 mA EW variation with rotational phase in AB 
Dor, and similar variations were recorded by Jeffries et al. 
(1994) for BD+22° 4409. These variations were interpreted 
in terms of rotational modulation caused by enhanced Li I 
from asymmetrically distributed starspots. It is tempting 
to speculate that the variability seen here is caused by the 
same phenomenon. However, it is then difflcult to under- 
stand (other than by appealing to the small sample size) 
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why no variation is seen on yearly timescales. It has been 
clearly demonstrated (for example by Barnes et al. 1998) 
that spot patterns on these young stars may only remain 
coherent for a month or less and in any case, it is unlikely 
that the stars were viewed at the same rotational phase (ro- 
tation periods are too uncertain to check this point) . A more 
likely explanation could lie in magnetic activity cycles. A 
cycle length of order 11 years (like the Sun) would mean 
that that any effects of plage or spot activity spread over 
the whole stellar surface may not show great variation from 
year to year, but possibly would over the course of many 
years or decades. There is evidence that some, but not all, 
cool stars share solar-like magnetic activity cycles (Baliunas 
et al. 1998). Few stars as active as the Pleiades G/K stars 
have been studied in such detail, although AB Dor's overall 
spot coverage does reveal changes on ~ 10 year timescales 
(Anders, Coates & Thompson 1992). 

A more sceptical interpretation is that much of the long 
term variation is caused by uncertainties in the setting of 
continuum levels prior to EW measurement. Even contin- 
uum level errors at the 1-2 percent level could result in mea- 
sured EW errors of about 10-20 mA. Effects of this size are 
certainly plausible given the S/N and different spectral reso- 
lutions of the datasets. Some support for this explanation is 
offered by the correlation between Ca I and Li I EW discrep- 
ancies in some, but not all, cases. It is therefore still possible 
that some of the variation is genuine. 

The results presented here also suggest that chromo- 
spheric activity variations on timescales of 1-10 years are 
limited to 20-30 percent in these active, young stars. This 
result is not entirely surprising, although some rare excep- 
tions were found by Soderblom et al. (1993b) in their more 
extensive Hq survey of the Pleiades. A similar lack of long 
term variability has been found in the coronal emission of 
Pleiades stars at X-ray wavelengths (Gagne, Caillault & 
Stauffer 1995), in contrast to the order of magnitude vari- 
ability in X-ray activity exhibited by the Sun over the course 
of its activity cycle. The lack of chromospheric activity vari- 
ability, combined with the order of magnitude dispersion in 
the chromospheric activity of the Pleiades stars considered 
here (presumably due to different rotation rates), make it 
difHcult to argue against active regions distributed over the 
entire surface being responsible for at least a proportion of 
the Li I EW scatter. Stuik et al. (1997) have modelled the Li I 
and Kl sensitivity to homogeneously distributed spots and 
plages on the surfaces of stars similar to those considered 
here. They find that Li I and K I EW scatters of ~ 100 mA 
are possible given the appropriate combinations of spot and 
plage filling factors. 

Stuik et al. (1997) have also shown that the Kl line is 
an extremely good proxy for the Li I line and mimics its re- 
sponse to different atmospheric stratifications. S93a noticed 
that there was a spread in Kl EW at a given B — V and 
this point has been stressed by Carlsson et al. (1994) and 
Stuik et al. (1997) as a reason to be sceptical about whether 
the spread in Li I EW genuinely reflects different Li abun- 
dances. The new data presented here conflrm this result and 
also that those stars with the largest Li I and K I EWs tend 
to be the most chromospherically active and rapidly rotat- 
ing. Indeed, S93a did not present Kl measurements for Hz 
174 and Hz 2407 which are among the most rapidly rotat- 
ing stars of their spectral type in the Pleiades, and have the 



highest K I EWs in this temperature range. These data rein- 
force the conclusion that there is a factor of two spread in the 
K I EW at a given B-V {in the range 0.8 < B-V < 0.95), 
compared with a factor of three spread in the Li I EW. While 
this Kl EW scatter remains unexplained, it is dangerous to 
attribute the spread in Li I EWs to abundance differences. 
Note though, that at cooler temperatures the spread in Li I 
EW grows while the Kl EW spread remains roughly con- 
stant, but the correlation between rotation and Li I EW also 
becomes much weaker. Perhaps this indicates, in agreement 
with the predictions of current models (Pinsonneault et al. 
1998), that during the PMS phase at least, rotationally de- 
pendent non-standard mixing is not the architect of spreads 
in the Li I EW. 



5 CONCLUSIONS 

I have observed a small sample of Pleiades late G and early 
K stars with consistent instrumentation and analysis tech- 
niques to look for variation in the strength of the Li I 6708A 
line on one year timescales. Variability on ten year timescales 
has also been searched for using previously published data. 
At the same time I have measured the strength of chromo- 
spheric activity and the K I 7699A line, which might plausi- 
bly have shown correlated variations. The following conclu- 
sions can be drawn. 

(i) The detection of Li I EW variability comparable to the 
spread in Li I EW in this temperature range would have im- 
mediately rejected the hypothesis that the spread reflected 
a genuine Li abundance scatter. I detect no variability on 
one year timescales that could correspond to to more than 
about 0.1 dex in the deduced Li abundances. 

(ii) On 10 year timescales there is some evidence for Li I 
EW variability which could produce 0.2-0.3 dex Li abun- 
dance spreads. Because these data were taken with differ- 
ent instruments and spectral resolutions, I remain sceptical 
about whether such variations are real, due to the difficulty 
in applying consistent continuum placements. 

(iii) The first two conclusions lead us to believe that what- 
ever causes the scatter in Li I EWs, it is not something that 
varies significantly on 1-10 year timescales. 

(iv) There is only 20-30 percent variability in the levels of 
chromospheric activity in this sample, despite a dispersion in 
chromospheric activity of an order of magnitude. This makes 
it difficult to rule out uniformly distributed magnetically 
active regions as a cause of Li I line strength scatter. 

(v) The data reinforce the conclusion that the strengths 
of both the Li I 6708A and Kl 7699A lines increase with 
rotation rate and chromospheric activity at a given colour. 
It seems unsafe to conclude that Li depletion in the effective 
temperature range 4800-5200 K is correlated with rotation 
until the large scatter in the K I EWs has been satisfactorily 
explained. 

A number of useful extensions to this work are apparent. 
Ideally, the sample size should be increased and extended to 
include cooler Pleiads. It would also be most advantageous 
to repeat the observations presented here in about 5 years 
time, using the same instrumentation if possible to avoid sys- 
tematic errors. Simultaneous photometry or doppler imag- 
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ing could be used to determine the actual surface coverage 
of spots or active regions at the time of observation. 
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